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^0 ' Patterned irradiation of cuprate superconductors with columnar defects allows a new generation 

, of experiments which can probe the properties of vortex liquids by forcing them to flow in con- 

' fined geometries. Such experiments can be used to distinguish experimentally between continuous 

disorder-driven glass transitions of vortex matter, such as the vortex glass or the Bose glass tran- 
sition, and nonequilibrium polymer- like glass transitions driven by interaction and entanglement. 
C ' For continuous glass transitions, an analysis of such experiments that combines an inhomogeneous 

scaling theory with the hydrodynamic description of viscous flow of vortex liquids can be used to 
infer the critical behavior. After generalizing vortex hydrodynamics to incorporate currents and 
^ ' field gradients both longitudinal and transverse to the applied field, the critical exponents for all 

six vortex liquid viscosities are obtained. In particular, the shear viscosity is predicted to diverge as 
^ ' \T — Tbg\~'^' at the Bose glass transition, with i/ ~ 1 and z ~ 4.6 the dynamical critical exponent. 

The scaling behavior of the ac resistivity is also derived. As concrete examples of flux flow in con- 
, fined geometries, flow in a channel and in the Corbino disk geometry are discussed in detail. Finally, 

the implications of scaling for the hydrodynamic description of transport in the dc flux transformer 
geometry are discussed. 
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I. INTRODUCTION 



In the mixed state of cuprate superconductors the magnetic field is concentratedja-an array of fiexible fiux bundles 
that, much like ordinary matter, can form crystalline, liquid and glassy phases.EJ u The dynamics of the flux-line 
, array determines the resistive properties of the material and has therefore been the focus of much theoretical and 
' experimental work. Of particular interest are the phase transitions connecting various forms of vortex matter. For 
. example, aJirst order melting transition is now believed to connect the crystalline Abrikosov flux lattice to a melted 
0^ flux liquid. El If the barriers to line crossings are sufficiently high, rapidly cooled vortex liquids can bypass the solid 
_ phase altogether and form an entangled polymer-like glass phase.a The transition to such a polymer glass would occur 
when the system gets trapped in a metastable state, much like the transition in ordinary window glass. 

Other types of glasses are also possible because of pinning in disordered samples. It was suggested some time 
ago that point disorder may dci^ie a continuous transition from a vortex liquid to an isotropic vortex glass state, 
with vanishing linear sesistivity.ElEl The existence of this transition is still in doubt, as many of the early experiments 
claiming to observe imu were in fact dominated by twin boundary pinning. The vortex glass remains, howevei, a 
credible candidate for the description of the disordered solid phase observed in untwinned crystals at high fields.EI 
\ Correlated disorder, that is disorder that can pin vortex lines coherently along a specific direction, is also very 
O important in many materials. It can be created artificially via the introduction of columnar damage tracks created 
by hea.vii_.ion irradiation - a procedure that was shown to result in dramatic improvement in the pinning of vortex 
lines. E3"E3 It can alsp-|be present in the material in the form of families of parallel twin planes that pass completely 
through the sample.Ej In samples with columnar disorder, if the concentration of damage tracks (assumed to pass 
completely through the sample) exceeds the number of flux lines, there is a low-temperature "Bose glass" phase, 
• in which every vortex is trapped on a columnar defect and the linear resistivity vanishes .Ej At high temperatures 
the vortices delocalize in an entangled flux-line liquid. The high temperature liquid transforms into a Bose glass 
via a second order phase transition at[3j»«. Samples with single families of twin boundaries also exhibit continuous 



transitions to anisotropic glass phases. 
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The Bose glass transition is the one that is best understood and experimentally characterized. For this reason we 
will use it in much of the following as our prototype of a disorder-driven continuous glass transition in vortex matter. 
It has been studied theoretically by viewing the vortex line trajectories as the world lines of two-dimensional quantum 
mechanical particlesoa. The thickness of the superconducting sample corresponds to the inverse temperature of the 
quantum particles. In thick samples the physics of vortex lines pinned by columnar defects becomes equivalent to 
the low temperature properties of two-dimensional bosons with point disorder. In the low temperature phase Bose 
glass phase the vortices behave like localized bosons. The entangled flux liquid phase is resistive and corresponds to 
a boson superfluid. For simplicity, we confine our attention here to fields less that the "matching field" = r7,pi,i(/>o, 
with ripin the areal density of columnar pins and (po = hc/2e the flux quantum. 

Although an exact theory of such continuous transitions from the Bose glass or from the vortex glass to the entangled 
flux liquid phase is not available, near the transition m ost jphjvsica l properties can be described via a scaling theory 
in terms of just two undetermined critical exponentaj'EiHiZHla^E^Eil. The low temperature disorder-dominated glass 
phase is characterized by a a correlation volume whose size diverges at the transition. Anisotropic disorder, such as 
columnar defects or twin planes, results in an anisotropic correlation volume. The most complicated case is that of 
planar disorder, embodied for instance by a single family of parallel twin planes. In this case there are three correlation 
lengths that diverge at the continuous glass transition temperature, Tq (see Fig. 1). The growth of correlations in 
the direction perpendicular to both the external field and the twin planes is described by the correlation length 

U{T)-\T-TGr. (1.1) 

A second diverging length, ^± (T) , describes the extent of correlations in the direction perpendicular to the external 
field, but parallel to the plane of the twin boundaries, 

|i(r)~ It-TgI-^ (1.2) 

Finally, there is a diverging correlation length along the external field (here the z direction), which is aligned with 
the twin planes, 

eil(T)^ei^|T-TGr''«, (1.3) 

with C the anisotropy exponent. For planar disorder = 1 + i>/iy (this constraint follows from the finiteness of cn at 
the transition). In materials with columnar pins aligned with the external field the correlation volume is isotropic in 
the plane normal to the columns. Thus = and there are two correlation lengths that diverge at the Bose glass 
transition transition temperature, Tbgj characterizing correlations in the directions perpendicular {S,±_) and parallel 
(5||) to the linear defects. In this case the anisotropy exponent ( has the value C = 2. The same results are expected 
to hold for a mosaic of twin boundary planes, all containing the field direction, on scales large compared to the 
mosaic size. Finally, for the case of point disorder coiitrolling the transition to the vortex glass phase at Tvg, the 
correlation volume scaling is assumed to be isotropic.El There is only a single diverging correlation length, £^±, with 
— 1, although it is hard to rule out different length scales parallel and perpendicular to the field direction. The 
correlation time controlling the relaxation of critical fluctuations is assumed to diverge as 

T^e^^lT-TG]-'", (1.4) 

with z the dynamical critical exponent. The values of the critical exponents v and z of course depend on the type 
of disorder and differ for instance for the Bose glass and the vortex glass. For the Bose glass transition the critical 
exponents have been determined via simulations to be ~ 1 and z ~ 4.6lj. 
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Figure 1. The two anisotropic disorder geometries discussed in the text and the corresponding diverging correlation 
lengths. 

Scaling can then be used to relate physical quantities to these diverging correlation lengths and time. In particular, 
the linear resistivity p±^ (T) of the vortex liquid for currents applied in the ab plane is predicted to vanish as T — > Tq 
from above as 

Pi(r)~ iT-rcr^^-?) (1.5) 



in bulk samples in three dimcnsioni 




Although some predictions of the scaling theory have been tested experimentally, there are as yet no direct mea- 
surements of the transport coefficients usually associated with glass transitions in conventional forms of matter, such 
as the shear viscosity. In this paper we show that the liquid viscosities also exhibit strong divergences at a continuous 
disorder-driven glass transition of vortex matter. We make explicit predictions for the critical exponents controlling 
the divergence of the liquid viscosities and propose experiments which test our predictions. For instance, the behavior 
of the shear viscosity is found to be determined by the dynamical critical exponent z that controls the divergence 
of the relaxation time in the glass phase. A direct measurement of the shmr viscosity would therefore provide a 
direct probe of the diwerging relaxation time associated with glassy behavior .E£l A brief account of this work has been 
presented elsewhere. c3 

As for ordinary ipatter, the shear rigidity of the vortex array can be probed by driving the vortices to flo w . in 
confined geometriesBJ The fabrication and use of such confined geometries was pioneered by Kes and coUaboratorsE^EHl 
to study the shear rigidity of the two-dimensional vortex liquid near the freezing transition in thin superconducting 
films. More recently, patterned irradiation of cuprate superconductors with heavy ions has made it possible to create 
samples with controlled distributions of damage tracks that will allow for a new generation of experiments that may in 
fact provide a direct probe of viscous critical behavior near the Bose glass transitioned. Specifically, by starting with 
a clean sample, at temperatures such that point disorder is negligible, it is possible to selectively irradiate regions of 
controlled geometry by covering the top of the sample with a suitable mask. Perhaps the simplest geometry that can 
be created is that of the channel shown in Fig. 2. This closely parallels an experimental setup proposed several years 
ago for measaring the shear viscosity of an entangled vortex liquid forced to flow in the channel between two parallel 
twin planes .Ej Now we propose that a similar geometry created by selective heavy ion irradiation may be used to 
probe the growing correlations at the Bose glass transitions 
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Figure 2. A weakly irradiated channel where the flux liquid is sandwiched between two heavily irradiated Bose- 
glass contacts. A current J applied across the channel yields flux motion along the channel. Vortex pinning at the 
boundaries propagates a length S into the channel, yielding a spatially inhomogeneous flow profile. 
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As shown in Fig. 2, the side regions of the samples have been heavily irradiated, and are characterized by a high 

(2) (2) 

matching field and transition curve Tj^^, while the channel is lightly irradiated with a lower matching field 

b'^^ < b'^^ and transition curve tI^I- When T^^^ < Tbg < T^bg' t^e flux array in the channel is in the liquid 
state, while the contacts arc in the Bose glass phase. Flow in the resistive flux liquid region is impeded by the "Bose- 
glass contacts" at the boundaries, as the many trapped vortices in these regions provide an essentially impenetrable 
barrier for the flowing vortices, due to their mutual interactions. As the temperature is lowered at constant field, so 
that the Bose glass transition T^q of the liquid region is approached from above (Fig. 3) the Bose glass correlation 
length increases, forcing the pinning at the boundaries to propagate into the liquid channel and yielding a spatially 
inhomogeneous electric field profile which can be probed experimentally. We will show that the spatial inhomogeneity 
of the electric field occurs precisely on length scales given by the correlation length, providing an unambiguous 
prescription for extracting critical exponents from this type of transport measurements. 

By controlling the disorder in the channel region, the same geometry can also be used to study the onset of 
rigidity in the vortex array near the vortex glass or the polymer glass transitions. In the former case the channel 
(now free p£ damage tracks) should contain a high concentration of point defects (augmented, perhaps, by proton 
irradiationEj) . The polymer glass transition may take place when the channel is "clean" , with only a low concentration 
of oxygen vacancies. A crucial difference between the polymer glass transition and the continuous disorder-driven 
glass transitions discussed above is that the former is a nonequilibrium phenomenon associated with the slowing 
down of the system dynamics on experimental time scales and does not exhibit universal critical behavior. For this 
reason, by testing for a divergent length scale, transport measurements in confined geometries may actually be usccL-tp 
distinguish experimentally between a continuous disorder-driven glass transition and the polymer glass transition.cIl 




Figure 3. A sketch of the {B, T) phase diagram for the fiux array in the weakly irradiated channel region. The heavy 
line, Bg^{T), denotes the continuous transition from the Bose glass to the entangled liquid. Also shown is the location 
B^g(T) of the Bose glass transition line in the heavily irradiated contacts. When a field B^qIT) < B < BgQ{T) 
is applied, the fiux array in the channel is in the liquid state, while the contacts are in the Bose glass phase. By 
decreasing the temperature at constants field along the direction of the arrow, the Bose glass transition of the channel 
region is approached from above. When B, instead of H, is used for the phase diagram, the Meissner phase collapses 
to a line, as shown. 

As shown in Ref. |2j, the electrodynamics of vortex liquids is in general nonlocal due to interactions and entan- 
glement, even away from phase transitions. Such nonlocalities may become very important near a phase transition, 
particularly when the transition is continuous. A natural way of incorporating the nonlocality in the long wavelength 
description of vortex liquid dynamics is via a set of hydrodynamic equations where the effect of intervortex interac- 
tions appears as viscous forces. A complete set of hydrodynamic equations for the vortex liquid is given in section ||. 
These equations incorporate driving forces from external current both parallel and transverse to the field H, as well 
as both compressional forces and pressure forces due to vortex segments that are not aligned with the z axis. They 
therefore generaliKe earlier hydrodynamic equations written down by us that only considered flow driven by currents 
in the ab planeEJ As pointed out earlier by Mou et alEa, the line nature of vortex matter yields to a proliferation 
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of viscosity coefficients, much like in liquid crystals. There are therefore six independent vortex liquid viscosities in 
these equations. The general equations are complicated and will be discussed in more detail below. The equation for 
the coape-grained flow velocity v(r) of the vortex liquid in the channel geometry of Fig. 2 is, however, very simple, 
namely£j 

- 7j_v + yyV^v -no0oz X J, (1-6) 
c 

where 7_l is the friction coefficient that incorporates both the interaction of the vortex cores with the underlying ionic 
lattice and the pinning by material defects and 77 is the vortex liquid shear viscosity. The term on the right hand side 
of Eq. ( |1.6| ) is the Lorentz current due to an external current J applied across the channel, with no the areal density 
of vortices. Since the local electric field from flux motion is e(r) = HoisLz x v(r), Eq. ( |1.6| ) can be rewritten as a 
nonlocal version of Ohm's law, 

e-(52vie = piJ, (1.7) 

where 

is the in-plane linear resistivity in a bulk sample. The viscous length 



5{T) - ^im/l^iT) (1.9) 

controls the importance of viscous drag. In a bulk sample and in the absence of strong inhomogeneities in the spatial 
distribution of disorder, the electrodynamic response is expected to be es sent ially local on large scales. Far from the 



boundaries in ma cros copic samples, the viscous drag drops out from Eq. (1/7) and the linear vortex liquid resistivity 



is given by Eq. (1.8) and is controlled entirely by the friction coefficient, 7_l. In contrast, in a conflned geometry 
viscous drag becomes dominant near a glass transition and the viscous length controls the flux liquid resistivity in 
this case. 

For continuous disorder-driven glass transitions, such as the vortex or Bose glass transitions, a dimensional analysis 
of the hydrodynamic equations, combined with the scaling hypothesis that the diverging correlation lengths and 
^11 are the only relevant length scales in the problem, leads to the identiflcation of the viscous length with the in-plane 
correlation length, 

5{T)^i^{T)^\T-TG\-r (1.10) 



Since the bulk resistivity vanishes at Tbg according to Eq. (1^), or equivalently the friction 7j^(T) diverges with the 
same exponent, 

l^{T)^\T^TBGr^''-^\ (1.11) 

the identiflcation of these two length scales immediately leads to the prediction that the liquid shear viscosity diverges 
at a continuous glass transition as 

ry(T) ^ ei(T)7±(T) ^ |T - TBGr'^*'+'"^^. (1-12) 

At the Bose glass transition, where C = 2 and ~ 1, the divergence of the shear viscosity is controlled by the 
dynamical exponent z. 

In contrast, there is-iuo obvious diverging correlation length controlling the polymer glass transition proposed some 
time ago by one of useI as an alternative to the vortex glass scenario and a possible explanation for the irreversibility 
line observed experimentally at high fields. In this scenario a viscous vortex liquid with large barriers to line crossing, 
when cooled rapidly may drop out of equilibrium at a polymer glass transition temperature Tpc, well before it has 
time to nucleate the equilibrium crystalline phase. The polymer glass is metastable and the transition is analogous to 
the glass transition conjectured in ordinary forms of matter. Disorder plays no role in the transition itself (although 
it may have dramatic effects in controlling the dynamical response of the vortex array) and the vortex friction 7j_ 
remains finite across the transition. The shear viscosity grows rapidly as T ^ ^pg temperature dependence 

might conceivably be described by the standard Vogel-Fulcher form, 

,7(T)«r;oe^/(^-^-), (1.13) 



5 



which works very well for glass transitions in conventional forms of matter. Experiments capable of measuring both the 
friction coefficient 7_l (which can be extracted from a bulk resistive measurements) and the viscous length S (which can 
be obtained from resistive measurements in confined geometries) can therefore distinguish between a kinetic polymer 
glass transition and a continuous-disorder driven glass transition, as well as providing direct measurements of the 
critical exponents in the latter case.EZI At a continuous glass transition, both friction and viscosity will diverge with 
a characteristic critical exponent. At a polymer glass transi tion, the friction will remain finite and the viscosity is 
expected to exhibit the Vogel-Fulcher behavior given by Eq. (1.13). 



As shown in section III, dimensional analysis of the general hydrodynamic equations yields predictions for the 



scaling exponents for all the six vortex liquid viscosities near a continuous glass transition. The results (see Sec. II 
for the precise definition of these quantities) are summarized in table 1. 
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Table 1. A summary of the critical exponents for the friction and viscosity coefficients (see Sec. II for definitions) 
near a continuous glass transition in three dimensions. Each transport coefficient diverges at Tbg ^ ~ where 
t =\T — Tbg\ITbg- The value of a for each transport coefhcient is indicated in the second column. The anisotropy 
exponent C should be chosen as Q — 2 for the Bose glass transition and C = 1 for the vortex glass transition. 



We stress that the identification of the dynamical length scale 5 with the static correlation length is by no means 
obvious. It mirrors the corresponding result obtained in the theory of continuous melting in two dimensions, where 
the liquid shear viscosity scales as 77 ~ VCd' with the average distance between, unbound dislocations, which in turn 
diverges at the transition, yielding a corresponding divergence of the viscosityE3 In this case the precise relationship 
between static and dynamical properties is one of the outcomes of the theory of two-dimensional melting which occurs 
via successive dislocation and disclination unbinding. For the vortex or Bose glass transitions the divergence of the 
correlation length is due to collective effects mediated by extrinsic quenched disorder. The assumption made here is 
that this growing correlations associated with disorder will also drive the divergence of dynamical properties in the 
liquid phase, namely the viscosities. 

This hypothesis has powerful implications. As shown in Ref. ^ and described in more detail in section |^ below, 
by combining scaling in spatially inhomogeneous geometries with hydrodynamics we are able to obtain not just the 
scaling exponents for the transport coefficients, but also the full form of the scaling function for the linear resistivity 
of the vortex liquid in a constrained geometry. The precise form of the scaling function depends on the details of the 
geometry and was given in ^ for the channel and for the Corbino disk geometry (Fig. 4) used recently by Lopez et 
al.Lj. Once the form of the scaling function is known, a quantitative analysis of transport experiments can in principle 
be carried out to extract the critical exponents. 

The plan of the paper is as follows. In section |l| we generalize the hydrodynamic equations of Ref. |2j to incorporate 
gradients and driving forces in both the directions parallel and transverse to the applied field. Upon neglecting the 
small transverse frictional forces responsible for the Hall effect, the hydrodynamic equations contain two friction 
coefficients and six viscosity coefficients. In section III, after briefly summarizing the scaling theory for the Bose glass 
transition in bulk, we discuss the dimensional analysis of the hydrodynamic equations and obtain the critical exponents 
for all the transport coefficients, as given in table 1. In section [V we combine scaling in spatially inhomogeneous 
systems with hydrodynamics to show how the full form of the scaling function can be obtained for specific geometries. 
The Corbino disk geometry is discussed in some detail. Finally, section V generalizes both the scaling theory and the 
hydrodynamics to finite frequency to discuss the critical scaling of the linear response of the vortex liquid to ac drives 
at the Bose glass transition. Appendix A makes contact with the earlier work by Mou et al.Ea on nonlocal effects in 
vortex liquids. In appendix B we apply scaling ideas to the analysis of transport in the dc flux transformer geometry. 
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II. HYDRODYNAMICS OF FLUX-LINE LIQUIDS 



Here we generalize the hydrodynamics of flux-hnc hquids discussed earher by us0 to incorporate the effect of driving 
currents in arbitrary directions, with components both parallel and perpendicular to the external field Hq = zHq 
responsible for the onset-af the vortex state. The resulting hydrodynamic equations contain the six viscosity coefficients 
discussed by Mou et alBa that control the nonlocal electrodynamic response of flux liquids. 

The electrodynamics of a type-II superconductor in the mixed state is described by Maxwell's equations for the 
local fields h{r,t) and e(r,t), 

-dth + Vxe = 0, (2.1) 

c 

Vxb=— j, (2.2) 

c 

where we have dropped the displacement current which is negligible at all but very high frequencies {lu ~ c/A ~ lO^^Hz) 
and j denotes the total current density, including the equilibrium response of the medium, in addition to any external 
current. In addition, the fields satisfy divergence equations, 

V-b = 0, (2.3) 
V • e = Anp. (2.4) 

Both fields can be obtained from a vector potential A and a scalar potential 0, according to b = V x A and 
e = —-dtA — V0. We work here with the Coulomb gauge, V • A = 0. The equations for the fields must be 
supplemented with a constitutive equation for the current j. In the linear (Ohmic) regime this is simply given by a 
nonlocal generalization of Ohm's law, 

Ja(r) =£sa/3(r,r')e/3(r'), (2.5) 

where Greek indices a, P, ... run over the values x, y, z and are used to denote Cartesian components in three dimensions 
and Sa/3 is a nonlocal conductivity tensor. Below we will also use Latin indices i, j, ... to denote Cartesian coordinates in 
the xy plane, i.e., Latin indices will only assume values x and y. Nonlocality must be incorporated when describing the 
mixed state of the cuprates where interactions and entanglement among the vortices can yield long-ranged correlations 



and large scale nonlocal electrodynamic response. It is useful for the following to rewrite Eq. (2.5) by sepa rating the 
current in an external and internal contributions, j = J + j'"', with j'"* = x (b — H). Equation (2.5) can then 
be written as 

J„(r)-c(Vx ^)^=£E„;3(r,r')e^(r'), (2.6) 

where — 47r|^ = b — H represents the local magnetization of the medium, with F{T,h) the Helmoltz free energy of 
the superconductor in an external field, H. The second term on the left hand side of Eq. ( |2.6| ) incorporates the static 
Meissner response of the material as well as the current due to local Magnus and pressure forces from the vortices. 

In this paper we are interested in fiuctuations in the b and e fields due to fiuctuations in the vortices degrees 
of freedom. The relationship between the local supercurrent, j^, and the vortex degrees of freedom is obtained by 
minimizing the Ginzburg-Landau free energy functional for fixed vortex configurations. Note that in general by using 
a two-fluid picture there will also be a "normal current" contribution to the internal current, j'"* — + in- We will 
not incorporate this normal part of the response here. In the London approximation, where the magnitude of the 
superconducting order parameter is assumed constant, and only fluctuations in its phase, 9, are retained, one obtains 
the usual London equation, 

^'j- = -|^(A-0oV0), (2.7) 

where A is the penetration length and — hc/2e the flux quantum. For clarity here we discuss only the case of 
an isotropic superconductor. The equations are easily generalized to a uniaxial material. In addition to the field 



fiuctuations described by Eq. (2.7), there are field fiuctuations representing thermal deviations from the solution of 
London equation, which are neglected here. 

In the chosen gauge (V • A = 0) the longitudinal part of the supercurrent is simply determined by the phase in the 
usual way, 
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(2.8) 



Vortices provide a source of vorticity for the supercurrent, as the curl of Eq. (2.7) gives 

c 



V X j. 



47rA2 



(2.9) 



where T is the vortex density vector, 



T{r,t) = —VxVe. 
Zn 



(2.10) 



We note that the curl of V6' is nonzero in the presence of vortices. By combining Eq. (2.9) with Eq. (2.2), we 
immediately obtain the London equation determining the field due to the vortices. 



b- A^V^b 



(2.11) 



To obtain the contribution to the electric field from the vortex degrees of freedom, we differentiate Eq. (2.7) with 
respect to time and use e = —-dtA — V^, with the result 



e = - V(/i H :^dtjs OtVe. 



The vortex part of the electric field is then given by 
where 



2c 



1 
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r(r,t) 



(2.12) 



(2.13) 



(2.14) 



The antisymmetric tensor Qg- , as d efined as in Eq. (2.14) describes the vortex current. It is precisely the vortex flux 
tensor as with the definition ( ^.14 ) the vortex density vector T satisfies the exact conservation law (conservation of 
vorticity) 



dtTa, + dpQap = 0. 



(2.15) 



The focus of the present paper is on the flux-line liquid regime of vortex matter. In order to take advantage of the 
vast phenomenology developed for describing long-wavelength static and dynamical properties of dense liquids with 
partcle-like or polymer-like degrees of freedom, we will focus our attention below on the vortex degrees of freedom, 
rather that on the fields b and e. When the z direction is chosen along the external field H, vortex configurations 
are conveniently parametrized in terms of a set of N single- valued functions rn{z), which specify the position of the 
n-th vortex line in the xy plane as it traverses the sample along the field direction. The three-dimensional position of 
each flux line is then R„(z) = [r„(z), z] and the vortex density vector and flux tensor are given by. 



n—1 

N 

Qapir, t) = Eap^ ^ [dzHn X v„^ 5{r - r„(z, i)). 



(2.16) 
(2.17) 



with r — (rj^,2:) and v„(2;,t) — (9fr„ the vortex velocity. The three-dimensional vortex density vector is also often 
written as 



(t,n). 



(2.18) 



where n is the local areal density of vortices and t describes the local tilt of vortex lines away from the z direction. 
Neglecting spatial inhomogeneities on scales < A, Eq. (2.11) gives 
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bz{r,t) = (j)on{r,t), 



Similarly, the vortex flux tensor Qa/s can be written as 

Qxx Qxy Qxz 
Qyx Qyy Qyz 
Qzx Q zy Qzz 



noV -nov^ 
—noV —riQVy 
hqVx noVy 



(2.19) 



(2.20) 



where v represent the flow velocity of vortices moving in the xy plane and uqV — ^^ijQij is the number of vortex 
segments aligned with the x direction crossing an area normal to the y direction per unit area and per unit time, with 



N 



nov(r,i) = ^ v„(5(r- r„(z,i)), 

Tl=l 

N 

naV{r,t) = (dzVn x v„ j(5(r - r„(z,t)). 



(2.21) 
(2.22) 



n=l 



We have explicitly linearized in the fluctuations as we are only interested in linearized hydrodynamic equations. 
Neglecting again nonlocalities on length scales ~ A, Eq. (2.13) yields 



2c 



(2.23) 



or 



(2.24) 



We are interested here in the long-wavelength properties of flux-line liquids on scales large compared to A and 
to the average intervortex spacing, oq = yj4>o/ Bq, with Bq the equilibrium mean value of the field, (b) = zBq- 
For this purpose we can abandon the description in terms of the microscopic vortex degrees of freedom 'R.niz) in 
terms of coarse-grained or hydrodynamic fields describing fluctuations in the conserved variables of the system and 
corresponding to those collective degrees of freedom whose relaxation rate vanishes in the long wavelength limit. For 
the vortex liquid the relevant conserved variables are the three components of the vortex density, T(r,t) = (t,n). 
The coarse-grained vortex density is defined as in Eq. (2.16), with the (5-function replaced by a smeared-out two- 
dimensional (5-function, 5bz{^a^), with a finite spatial extent of the order of the inverse of the Brillouin zone boundary, 
ksz = \/47mo, 



3SZ 



(2.25) 



q±<kE 



As discussed elsewhere, the long wavelength equilibrium properties of vortex liquids can be described in terms of a 
coarse-grained free energy that fully incorporates the nonlocality of theiiitervortex interaction. Such a free energy can 
be derived by explicit coarse-graining of the microscopic vortex energycil or can be written down phenomenologically 
using familiar ideas from liquid state physics. It takes the form of an expansion in the fluctuations of the hydrody- 
namic fields from their equilibrium values and of the corresponding gradients. Retaining only terms quadratic in the 
fluctuations, it is given by 



2nl 



CL{r - r')(5n(r)(5n(r') + C44(r - r')t(r) • t(r') 



(2.26) 



no and cl and C44 are the nonlocal compressional and tilt elastic moduli of the flux liquid, given 



where 5n ^ 

elsewhere. Eil 

The vortex density vector, T — (t,n), satisfies an exact conservation law, Eq. ( 2.15| ), or, in terms of n and t, 



dfU + V • ngv = 0, 

dtU + djQij = dzUoVt, 



(2.27) 
(2.28) 
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with the condition V • T = 0, or 



d^n + V_L-t = 0, (2.29) 

ensuring that flux lines do not start nor stop inside the sample. 

To obtain a closed set of hydrodynamic equations it is necessary to specify the constitutive equations that express 
the flux tensor Qa/s in terms of the density field and its gradients. Since the vortex dynamics is overdamped, the 
constitutive equation is simply ob tain ed requiring that the net force on each local volume of flux liquid vanishes. Its 



form is analogous to that of Eq. (2.6) and it is given by 



Itiv.XpQxp — Vfiu,af3,XpdadpQxp + "^j^ ^ ^t" ^ ~^pi^xJx- (2.30) 

The first term on the left hand side of Eq. ( ^.3C| ) represents the frictional force per unit length on the vortex liquid, 
with 7;ji/,Ap a friction tensor. The second term describes viscous forces from intervortex interaction and entanglement, 
with f}ij_u^a(3,Xp the viscosity tensor. The third term contains Magnus and pressure forces arising from vortex density 
gradients and the last one is simply the Lorentz force due to an external current, J. 

Neglecting for simplicity transverse drag forces associated with the Hall effect, there are two independent components 
of the friction tensor and six viscosity coefficients. It is useful to rewrite the constitutive equations in terms of the 
in-plane vortex flow velocity v and the tilt current V , 



7-LV = [r/Vi + 7hdl]^+ribV ^{V i_ ■ v) - 77^ (z9, x V i_)V 

SFl g SFj 
5n ^ St 



"■0 V±— d,^—] zxJ. (2.31) 



711 = [r7*Vi + vld^,]V^ ri^id, ■ (Vj. x v) + noz 



5Fl 
5t 



nQ(j)Q 



J. 



(2.32) 



The constitutive equations (2.31) and (2.32) are now written in terms of forces per unit volume. The coefficients "f±_ 



and 7|| have dimension of friction per unit volume, as usual, and 77, 77^, ?7f,, 77 , 77* and 77^; have dimensions of fluid 
viscosities. 

The six viscosity coefficients appeamig in Eqs. (2.31) and ( 2.32| ) are related to the six components of th e nonloca. 1 
cond uctivity d iscussed by Mou et alE3 The precise relationship is displayed in Appendix A. Equations ( 2.27 - 2.29 ) 
and (2.31,2.32) provide a closed set of hydrodynamic equations to describe the response of vortex liquids .to external 
currents and fields. These equations generalize equations written down earlier by uscj and by othera£3'E2l. by incor- 
porating all pressure forces, including those from vortex segments that are not parallel to the external field direction, 
all nonlocal viscous forces, as well as external forces due to driving currents aligned with the external field. In the 
next section we will show how predictions for the singular behavior of for the six viscosities near a continuous phase 
transition from the vortex liquid to a glassy state can be obtained by simple scaling arguments. 



III. SCALING THEORY 



The divergence of pinning energy barriers for vanishing driving currents underlies both the theories of the Bose 
glass and the vortex glass transitions, which are expected to be continuous. The properties near the transition have 
therefore been described by a scaling theory in terms of divergent lengths and time scales. A scaling theory for the 
transition from a vortex liquid to an isotropic vortex glass was developed by Eisher, Fisher and Husdj and then 
adapted by Nelson and Vinokur to the transition to the anisotropic Bose glassJlJ Here we summarize this scaling 
theory in a unified way that applies to both transitions by keeping track of the anisotropy exponent C that has value 
C = 2 for the Bose glass and C, = 1 for the vortex glass. Of course the values of the other critical exponents are also 
different in the two cases. 

The (Gibbs) free energy density g of the vortex array is a function of temperature T and external field, H — 
zi^ii + H^, with g = g{T, H±, )• Here = Ha + 6H^^ , where Hq is the field responsible for setting up the vortex 
state and the flux lines are on average aligned with the z direction. As predicted by the renormalization group, 
the singular part of the free energy density, g^, is assumed to obey a homogeneity relation at the transition where 
temperature and the external fields are rescaled by different factors, 

gS,H^,5H^) = l-^''-^+<^gs{l^/''t,l^^H^,l^i^SH^^), (3.1) 
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for arbitrary length rescaling parameter, I. Here t is the reduced temperature t = \T — Tg\/Tg, with Tq either the 
Bose glass or the vortex glass transition temperature, and d denotes the system dimensionality. The case of interest 
here is d = 3, corresponding to a bulk superconductor. In this section A_i_ and Ay denote critical exponents, not to 



be confused with penetration lengths. Making the usual choice / = t^'^ 
relate physical quantities to the diverging correlation lengths and ^|| 
The local field in the superconductor is given by 



S^±, dimensional analysis will allows us to 
as described in Ref. |jl4|jl^Jl||^,|l|. 



-4n 



59_ 
6U 



and therefore scales according to 



'II 



^ ,-(d-l+C+Aj_) 
;-(d-l+C+A||)_ 



The longitudinal and transverse components of the response function are 

^^9s , nl(t)l 



+ 47rx±) = 



1 



47rx|| 



sm cii ' 



C44 



(3.2) 



(3.3) 
(3.4) 



(3.5) 
(3.6) 



where xi- ^-nd x\\ ^-re the components of the magnetic susceptibility and C44 and cn the tilt and compressional elastic 
moduli of the vortex array. The scaling of the susceptibilities is then given by 



XII 



C44 
1 

Cll 



^ ^-(d-l+C+2Ai) 
p(<i^l+C+2A||) 



(3.7) 

(3.8) 



At the vortex glass transition the response of the vortex array is rotationally invariant and both susceptibilities are 
expected to remain finite. This requires 



A!^ = Ap = -(d- 1 + 0/2. 
The scaling of the local fields at the vortex glass transition is then 

1 



Jd-i)/2A/-i 



Zv/2 



(3.9) 



(3.10) 



where the last equality applies for d — 2i. 

At the Bose glass transition, in contrast, the tilt modulus is expected to diverge as C44 
finite. In this case the response of the system is anisotropic, with 



\BG 



-(d-l)/2, 
-(d- 1 + 0/2, 



Z^, while Cll remains 

(3.11) 
(3.12) 



and the fields scale as 



Jd-l)/2A/2 



(3.13) 



where again the last equality in each of Eqs. (B.13) applies for d = 3. 

To determine the scaling of fields and currents that control the resistive properties of the vortex array, we note that 
gauge invariance of the Ginzburg-Landau theory requires that the fluctuating vector potential scales according to 
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A, 



The scaling of currents and fields is immediately obtained from their definitions, J — df /dA and E = 
where we denote by E the macroscopic (spatially-averaged) field from vortex motion, ft is given by, 



Ji 



and 



El 

En 



(3.14) 

-(l/c)dA/dt, 

(3.15) 
(3.16) 



A scaling ansatz for the IV characteristic can be written down based on the structure of the RG flows as 

Z(<i-2+C)j, 



E^{T,Ji_)^l-^^+''^E^(l^/''t,- 

i?ii(r,j||) = r(c+^)£;,i(/v^t,^ 



cknT 



With the choice I = t " ^ C^{T), we obtain 



E±ihJ±)^i 



d-2/ 



' cksT 



(3.17) 
(3.18) 

(3.19) 
(3.20) 



where £^ || are scaling functions. From these one can obtain various results described in the literature. In particular, 
by linearizing Eqs. (3.17) and ( [3.18 ) in the ohmic regime, one finds that both components of the linear flux liquid 
resistivity vanish as the transition is approached from the liquid side, according to. 



P||(T) = 



no 00 



no 00 



1 



'■ 1 



IT-Ti 



BG 



(3.21) 
(3.22) 



Equivalently, Eqs. ( 3.21 ) and ( [3.22| ) give the scaling exponents controlling the divergence of the longitudinal and 
transverse friction coefficients, 7||(r) and 7^(T). 

Dimensional analysis of the hydrodynamic equations now allows us to extract predictions for the scaling of the 
six viscosity coefficients at the transition. The scaling dimension of the fluctuating density and tilt fields n and t is 
t he sa me as that of the components of the local field b and therefore is obtained immediately from Eqs. ( 3.10 ) and 
(3.13) for the vortex and Bose glass, respectively. We scale all length and time scales in the hydrodynamic equations 
with the correlation lengths ^± and ^|| , using the scaling dimension just discussed for the various physical quantities 
appearing in these equations. By requiring the equations to be scale invariant near the transition, we can then read 
off the critical scaling of the viscosity coefficients. 



rj 

Vb ' 

Vz 

rf. 

Vx 



- ei7ii 



z+5-d 



i 

^ s_L 



^z-\-3 — d^ 



(3.23) 



The critical exponents for all the friction and viscosity coefficients are summarized in Table 1. 
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IV. SCALING AND HYDRODYNAMICS IN CONFINED GEOMETRIES 



The flow of vortex liquids in confined geometries provides a powerful experimental tool to probe the shear resistance 
of voEteps-jnatter. Experiments of this type have been carried out by the group of Peter Kes in thin superconducting 
films.E3^c3 In this two-dimensional case the vortex array is expected to melt from a solid into a liquid via the dislocation 
unbinding mechanism.E3 The shear modulus of the lattice vanishes continuously at the transition and is replaced by 
a shear viscosity in the liquid phase. The shear viscosity diverges as the transition is approached from the liquid 
side. In order to probe this behavior, Theunissen et al. studied the flow of a two-dimensional vortex solid through 
narrow channels in a double-layer device, consisting of a weakly pinning amorphous Nb^Ge bottom layer, covered 
with a strongly pinning NhN thick layer. Using nanolithography, parallel channels were etched through the NbN 
layer into the NbGe layer. By applying a current normal to the direction of the channels, vortices where-then forced 
to flow along the weakly pinning channels, while the NbN regions provide the strong-pinning boundaries.E3 By fitting 
their data to the form-for the flux-flow resistivity of a vortex liquid in a channel obtained earlier by us via a simple 
hydrodynamic modeO, these authors were able to extract the temperature dependence of the shear viscosity near 
melting, finding a behavior in agreement with the predictions of the theory of two-dimensional melting.cj 

As discussed in the Introduction, patterned irradiation of cuprate |-superconductors with heavy ions allows for 
experiments similar in spirit to the one carried out by Theunissen et al.B3. These experiments can probe the growing 
shear viscosity of the vortex array near a continuous glass transition and may shed some light on the connection 
between static and dynamical properties of these glassy systems. 

Scaling theory can again be used to determine the temperature dependence of fields and currents in this spatially 
inhomogeneous situation arising from the confined geometry. A generalized homogeneity relation for the local electric 
field e±_ perpendicular to any correlated disorder from flux motion at position a; in a channel of thickness L takes the 
form 

e,(T, J„ L) = r (^-)e, '^^^^^^^ f f ) ^ (4-1) 

By choosing again I = t~'^^ ^ we obtain 



In the entangled flux liquid the response is linear at small current. Upon expanding the right hand side of Eq. (4.2 
we obtain for Jj^ ^ 



e^{J^^Q,x,L)c^ piles' ^J^T{x/U,L/U), (4.3) 



where p*]^ — ^riQ^o/c^ (I/Ti,) ^^i^ Bardeen-Stephen resistivity of noninteracting flux lines, with 7j[ a bare friction. 

A scaling form for the resistivity p±{T,L) = AV/ {LJ±), with AV the net voltage drop across the channel, is easily 
obtained by integrating Eq. (p^), with the result. 

p^{T,L)=p^{T)f{L/U) (4.4) 

/2 



with f{x) = i J^l^ duJ^{u, x) a scaling function and p±{T) the bulk flux liquid resistivity given in Eq. (1^). Here j± is 



the renormalized friction coefficient that incorporates the growing Bose glass correlations near Tbg, 7_l = 7 m — ) , 

\ °o / 
,^ , ■ For L '3> the channel geometry has no effect and one 

must recover the bulk result. This requires f{x^l)^l. 

The scaling function can be determined by assuming that the long wavelength electric field of Eq. (|4.3| ) is 
described by the set of hydrodynamic equations discussed in Ref. For simple geometries where the current is 
applied in the ab plane and the flow is spatially homogeneous in the z direction, these reduce to the single equation 



for the coarse-grained flux liquid flow velocity v(r) given in Eq. (1.6). Intervortex interaction at the Bose-glass 
boundaries yields "adhesion" of the flux liquid to these boundaries, forcing the flux flow velocity to vanish or become 
very small. In the hydrodynamic model this translates into a no-slip boundary condition for the flux liquid velocity 
at the boundaries. By preventing the free flow of flux liquid, the Bose glass boundaries can significantly decrease the 



macroscopic fiux-fiow resistivity of the superconductor. The electric field profile is obtained by solving Eq. (1.7) for 
the appropriate geometry and with suitable boundary conditions. 
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A. Channel geometry 



The solution of the hydrodynamic equation for the simple channel ge omet ry sketched in Fig. 2, with a homogeneous 



current J — —xJ± apphed across the channel has been given elsewhere 
is given by 



e±{x,L) = pj_Jj_ 



1 



cosh(a;/^^) 
cosh(L/2e±) 



and is displayed here for completeness. It 



(4.5) 



Upon comparing Eq. ( |4.5| ) to Eq. (4.3), we identify the quantity in square b rack ets in Eq. (4.5) with the scaling 
function J^. The scaling form for the resistivity is obtained by integrating Eq. (4^), with the result 



1 



2^ 
L 



tanh 



L 



If <C we recover the bulk result given by Eq. ( 3.21 ), 



(4.6) 



(4.7) 



Near the transition, however, ^ L, and the resistivity depends on the channel width and is controlled by the shear 
viscosity, with 



P±l{T,L) 



(4.8) 



The strong divergence of the viscosity implicit in Eq. (48)is precisely the kind of behavior expected at a liquid-glass 
transition. The Bose glass transition is an example of a glass transition that is well understood theoretically and 
where precise predictions are available. 



B. Corbino disk 



Experiments in patterned geometries near the Bose glass transformation provide an exciting opportunity to probe 
viscous behavior near a second order glass transition. Other patterned geometries can be designed and used to probe 
some or all of the other viscosity coefficients. Of particular interest is the Corbino disk geometry, recently used by 
Lopez et al. for defect-free materialsE3 to reduce boundary effects in the flux flow measurements. This geometry was 
briefly discussed in ^ and is sketched in Fig. 4. 




Figure 4. Top view of the Corbino disk geometry with Bose glass contacts. The magnetic field is out of the page. 
The vortex array is in the Bose glass state in the inner and outer densely dotted regions and in the flux liquid state 
in the weakly irradiated annular region. A radial driving current drives flux motion in the azimuthal direction. 
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In ^ we proposed the fabrication of a Corbino disk with Bose glass inner and outer contacts. The annular region 
corresponding to Ri < r < R2 is in the vortex liquid state. A constant current / is injected at the inner boundary of 
the cylindrical flux liquid region and extracted at the outer boundary. The resulting current density is radial, with 



/ f 



(4.9) 



with W = R2 — Ri the width of the annular flux liquid channel. This current induces tangential flow of the vortices, 
which in turn yields a radial voltage drop that can be probed by a suitable set of voltage taps, placed as sketched in 
Fig. 4. If the intervortex spacing ag in the flux liquid is small compared to the size R2 — Ri of the annular channel 
where flux flow occurs, the flow can be described by the hydrodynamic equation (1.7). Using cylindrical coordinates 
with the z axis directed along the direction of the external field H || c, the vortex flow velocity is in the azimuthal 
direction, v — Vc/,{r)4>, corresponding to a radial field, er{r) — ^^^^^v^{r). The nonlocal Ohm's law takes the form 



r clr ar \ / 



(4.10) 



where ul = no(l)oI / {2ttWc) controls the strength of the Lorentz force. The general solution of Eq. 4.1C is given by 



Crir) 



P±I 



2TrW^A 



6 



(4.11) 



where Ii and Ki are Bessel functions. The constants ci and C2 are determined by requiring the liquid flow velocity, 
and therefore the electric field, to vanish at the boundaries with the Bose glass regions, er{Ri) = er{R2) = 0, with 
the result. 



Cl 



C2 



^^^1(^2)^1 - Ki{pi)/p2 



K,{p,)h{p2) 

h{pi)/p2 



Kl{p2)h{pi) 

h{p2)lpi 



(4.12) 



KMh{p2) - K^{p2)h{piy 



where pi 2 — Ri,2/S.±- As in the channel problem, spatial inhomogeneities are controlled by the in-plane correlation 
length ^j^. The field profile for a few values of in a micron-scale sample is shown in Fig. 5. 




Figure 5. The electric field profile for the Corbino disk geometry, with i?i = 2pm and R2 = \2pm, with an annular 
width R2 — Ri = 10pm. The vertical axis represents the reduced field 2T:er{r)/{p±I) in pm. The dashed line is the 
electric field for vanishing viscosity, e°(r) = 27T{Rt-Ri) r • ^^'^ shaded regions represent the Bose glass contacts. 

The field nonlocality arising from the growing correlation length near the Bose glass transition can be probed 
experimentally by measuring the voltage profile across the annular channel. A series of equally spaced contacts 
are placed radially, as shown in Fig. 5, at positions r„ = i?i + nd, for n — l,2,...,N. For the magnetic fields 
and temperatures of interest, the spacing d of the contacts is large compared to the intervortex spacing oq, so that 
hydrodynamics can be applied. The voltage drop across two successive contacts is obtained immediately by integrating 



Eq. (4.11) and it is given by 
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Ay(r„+i,r„)= / drer{r 



(4.13) 



2ttW 



I ln(r„+i/r„) + ci [/o(r„+i/^^) - /o(r„/C±)] 



-C2[ifo(r„+i/e±)-/^o(r„/ei)]}. 



The temperature dependence of the correlation length can then in principle be extracted by fitting the measured 
voltage to Eq. ( [1.14| ) . 

For comparison it is useful to discuss the limiting cases of very small and very large correlation length. Away 
from the Bose glass transition, where << -Ri,i?2,d, viscous forces are negligible and the boundary condition of 
vanishing flow velocity at the Bose-glass boundaries cannot be satisfied. The field from vortex motion is simply 



P±I 1 
2TrW r' 

The voltage drop across each pair of contacts is then given by 



AVb(r„+i,r„) 



2nW 



ln(r„+i/r„) 



(4.14) 



(4.15) 



and is controlled by the bulk flux flow resistivity p± given in Eq. (1.8). In this limit the vortex liquid only has very 
short ranged correlations; the local electric field is simply proportional to the local current and decreases as ~ 1/r. 
As a result, the voltage across every successive pair of contacts decreases monotonically at all temperatures as one 
moves from the center to the outer perimeter of the Corbino disk, i.e., Vo(r'2,ri) > Vo(r'3,r2) > ... > Vb(r'„+i, r„). 

Near the transition, when the correlation length becomes comparable to the channel width, viscous effects are 
dominant and the field profile can be approximated by 



with 6i, ^2 > given by 



e-(r) 

bi = 
b2 = 



P±I 



bir r ln(r/^_L' 

r 



Rl HR2/U)-Rl ln(i?i/6 



d2 7j2 



7?2 I?2 

U2 - Hi 

2 ln(i?2/i?i) 



(4.16) 



(4.17) 



The voltage drop across two successive contacts is then given by 

P±I f 26i 



A14o(r„+i,r„) 



y n+l '11 J 



bo In (r. 



i/rn 



1 



In (r„+i/^±) - rl In (rn/U) 



(4.18) 



The first term on the right hand side of Eq. 4.16 can be easily understood as it corresponds to the field resulting from 
the rigid body rotation of^a. vortex solid at an angular frequency lu = no(poI / {SircWr]). In this rigid body limit the 



electric field grows with rH3 As a result, the voltage across two successive contacts increases at all temperatures as 
one moves from the center to the outer perimeter of the Corbino disk, i.e., Vo{r2,ri) < Vo{r3,r2) < ... < Vo(rn+i,rn). 

In general as the temperature is decreased towards the Bose glass transition, one will observe a crossover from 
liquid-like (e^ ~ l/?") to solid-like (e^ ~ r) response of the vortex array. 

The voltages across each successive pair of five contacts are plotted in Fig. 6 as functions of d/£,±{T). The correlation 
length £,±{T) increases with temperature as the Bose glass transition is approached from above. The figure therefore 
displays the g rowth of the voltage with increasing temperature. At high temperature, where £,± << d, the voltage is 
given by Eq. ( [4.15 ), which is scaled out, and all curves approach 1. The scaled voltage drop across successive pair of 
contacts is no longer a monotonic function of temperature. Rather it first increases as then decreases as one moves 
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from the center towards the outer border of the disk. 



1.1 
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ll 
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Figure 6. The voltage drop across a series of equally spaced contacts located radially at r„ = Ri+nd, for n=l,2,...,5. 
The vertical axis represents the voltage across two neighboring contacts scaled by the corresponding voltage obtained 
for ^± « d, i.e., K+i,n = Al/(r„+i,r„)/AVb(r„+i,r„), where Ay(r„+i,r„) and AVb(r„+i,r„) are given by Eqs. 
(4.13) and (4.15), respectively. We have chosen Ri = d and R2 = Qd. The voltage is plotted as a function of d/^± 
which increases with temperature, as discussed in the text. 



V. RESPONSE AT FINITE FREQUENCY 

The scaling theory of the continuous Bose glass transition in bulk samples is easily generalized to a finite frequency 
, as described in Ref. [l^. This is simply done by the addition of another scaling combination ~ to the scaling 



fimctions of Eqs. (3.17) and (3.18). Similarly, we can generalize to finite frequencies the scaling theory for spatially 
inhomogeneous flow in constrained geometries. For simplicity we refer again to the channel geometry of Fig. 2 and 
only discuss the response of the vortex liquid to a spatially homogeneous current of finite frequency cu applied across 
the channel in the ab plane. The generalized scaling ansatz based on a renormalization group flow for the spatially 
inhomogeneous field from flux motion is. 



(T, Ji, X, L) = r(i+^)ei(;i/'^i, jj^, PcjTo, r^x, r^L), (5.1) 



with To a microscopic time scale.L3 The flow diagram is sketched in Fig. 7. Clearly the behavior near the critical 
point depends on how the transition is approached in this plane. For a; = we flx the length I as before by choosing 
^ til''' = 1, which gives ^±(r) ^ = t ^ and yields the temperature dependence discussed earlier for the dc 
response near the transition. The same analysis should hold at finite, but low frequencies, provided at we have 
t(Z*) >> with u)(Z) = ujtqV . This corresponds to the condition u) « \/t{T), with t{T) ~ the relaxation 

time of critical fluctuations. 
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Figure 7. A sketch of the flow diagram in the (w, t) plane for a few initial values of ajTo/f^^. 

Conversely, we can approach the transition along the frequency axis by letting T — Tbg and using a finite frequency 
to tune the distance from the critical point. In this case we fix I by choosing = 1, which gives ~ = w"^/^. 
If the transition is approached along the T = Tbg line, the scaling ansatz becomes, 

e±(T = Tbg. ^, x, L) = {^y^'+^^E^ ({^17'^'+^^ J u {Ci_)-'x^ it^r'L) . (5.2) 

As before, we are interested in the response in the vortex liquid phase, where the current-voltage characteristic is 
linear at small currents and we can write 

e^{T^TBG,J±^0,u;,x,L) = p"^(^^^ J^T^(x/CL,L/CL). (5.3) 

By integrating the field across the channel we immediately obtain the net voltage drop and the corresponding ac 
resistivity as, 

P±l{Tbg,L,uj) = p_l(Tbg,w)^ {l/Ca), (5.4) 

where 

P^{Tbg, ^) = pI - (^c.)l-2/^ (5.5) 

is the ac resistivity in bulk samples.0 Both real and imaginary parts of the linear frequency-dependent resistivity 
vanish as w — > with the same critical exponent, as obtained earlier in Ref. The scaling function in Eq. ( |5.3| ) 

is simply fui{x) = ^ J^^j^ duJ-'i^{u,x). For x >> 1 the ac resistivity in the channel must reduce to the bulk result, 
which requires fu^ix >> 1) = 1. 

To determine the form of the scaling function we assume again that at long wavelengths the local field from flux 
motion is described by the hydrodynamic equations discussed in section O. T he use of hydrodynamics for describing 
the ac response of current-driven vortex liquids has been discussed in Ref. |35|. Both friction and viscosity coefficients 
are in general frequency-dependent. For the simple channel geometry, there is again only one hydrodynamic equation 
for the Fourier components of the the flow velocity at frequency w, given by 
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7_l('^)v + rj{ijjy\/ V — -no(/)oz x J(^)- 
c 



(5.6) 



Determining the frequency dependence of the friction and viscosity coefficients requires a microscopic calculation 
that is beyond the scope of the present pap er. However, their scaling with frequency at the transition can be readily 
obtained. The hydrodynamic equation ( ^.6[ ) defines a frequency dependent viscous length, ^ ri(uj) / •f ±{uj) . It is natural 
to identify this length with the diverging Bose glass correlation length, 



(5.7) 



The complete scaling function defined by Eq. (5.3) is then readily obtained by solving Eq. (5.6) for no slip boundary 
conditions and it is found to have the same functional form as that obtained for the dc response. The ac resistivity 
for the channel geometry if given by 



Pj_l{Tbg,L,uj) = /9_l(Tbg,^) 



1 - ^ tanh — — 



(5.8) 

where p±{Tbg,'-^) is the bulk ac resistivity given in Eq. ( ^.5[) . As for the dc case, one obtains the following behavior. 



\ c J 7_L(a;) 
P±l(Tbg,L,uj) ['~^) 



12ri{uj) 



{iujy-'^\ for « L, 

^ L^iilY^ - L^(%^), for i'l«L. 



(5.9) 
(5.10) 



The condition ^ L defines a sample-dependent crossover frequency iol = {o-o/LY /to ~ below which the 
channel geometry becomes important and p±L ~ itoL'^. These results describe a Meissner-type response of the 
superconductor that arises because the vortices are pinned and immobile in the Bose glass phase. 

The behavior observed as the critical point is approached along an arbitrary line in the (w, t) plane depends on 
the location of that line. In general we expect that the ui = scaling results apply for ui << l/r(T) ^ t^, while the 
scaling for T = Tbg will apply for w >> l/r. The various regions for the expected scaling of the ac response are 
sketched in Fig. 8. Here the temperature scale is defined by £,±{T) ~ L, with ^ {ao/Ly/". For t > ti, the vortex 
array is insensitive to the confined geometry and the response is that of a bulk sample. For temperatures below tL, 
the geometry becomes important. 
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Figure 8. A sketch of the scaling behavior expected for the ac response in various regions of the {oJ,t) plane. The 
solid line is a plot of uj — 1/t{T) = t^" /tq. The frequency and temperature scales u)l and represent the scales 
below which the confined sample geometry becomes important. They are defined by = L, or ujltq = {ao/L)^, and 
by £j_(Ti) = L, which gives ~ (oq/L)^/'^. The ac resistivity p±l(uj,T) is expected to scale as follows deep in each 
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of the four regions sketched in the figure: region I, p^L ~ L'^t'^^; region II, pj^L region III, p^L {i'^y 

region IV, 
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APPENDIX A: COMPARISON WITH THE HYDRODYNAMIC FORMULATION OF MOU ET AL. 

Mou et al. have used time-dependent Ginzburg-Landau theory to evaluate the nonlocal dc conductivity tensor 
S^^(r,r') defined in Eq. (2.5). Assuming nonlocal effects are slowly varying on the length scale of interest, one can 
expand the right hand side of Eq. (2^) retaining terms up to quadratic in the gradients of the electric field. At long 
wavelengths, the nonlocal relation between fields and currents then takes the form. 



j^(r) = (T^^(O)e^(r) - S^api^dadpe^ir), 



(Al) 



where a^,y{0) is the usual long- wavelength conductivity and the tensor Sp,a(3u describes the effect of nonlocality to 
quadratic order in the gradients. This fourth order tensor has six independent components that were evaluated in 
Ref. ||. 

Upon comparing Eq. (A.1) to Eq. ( 2.30| ) and recalli ng th e relat ionsh ip between the fields and the vortex densities 
and fluxes given by the long- wavelength limit of Eqs. (2.11) and (2.13), 



b = 0oT, 



(A2) 



we see immediately that the local conductivity tensor cr^^(O) and the tensor S^apv are simply proportional to the 
friction and viscosity tensor, respectively, according to. 



00 



00 



(A3) 
(A4) 



Neglecting Hall voltages, the local conductivity tensor is diagonal and has only two independent components, due 
to the rotational symmetry in the plane normal to the applied field {xy plane), (TxxiO) — ayy{0) and (Tzz(O), which 
gives 



where S^^ 



cra/3(0) = axx{0)S^p + a^,{0)SazSf3^, (A5) 
Sa±,i3± = Sap — SazSf3z and the labels with a _L subscript, e.g., a±,f3±, only run over the two values 



X and y. The friction tensor is given by 

_ 1 / 



(A6) 



with 



CTxxiO) = 



nQ(f)a 



^..(0) = (^y 



7±, 



7||- 



(A7) 



Symmetry under rotations about the z direction also restricts the number of independent components of the viscosity 
tensor and the related tensor Sf^a^u to six. The tensor S^j,ai3v is given by 
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Sf^afSv — Sxxxx5'^ij5al3 + {Sxyyx Sxxxx) S'^^^ap ~ 2^^^^"^P'' ^ ^pP'^ai') 
+ SxzzxSjj"'6^^ + SzxxzS^i^S^^'S''^ + SzzzzSpzSczSpzS.z 
+ Sxxzz{SazSfj,z5p^ + 5i3zS,yzSap + 5/3zSfj,zS^„ + 5azS„zSpfj) ■ 

The viscosity tensor is given by 



with 



2nl 



+VzSazSpz^pu'rj_£'Yj_\p + 11^ap^tJ.vz<^z\p + V z^ az^ [3z^pvz^z\p 



o _ ( c \^ 

^xzzx — 1 , J ^Z 



{Sxyyx Sxxxx) — ( ] ) ^67 



Sz 
Sz 



nQ(j>o' 



^xxzz — n\ I ] ^x- 

2 V 71000/ 



(A8) 



(A9) 



(AlO) 



Mou et al. evaluated the lowest order fluctuation contribution to the components of the tensor S^apv for a vortex 
liquid using a time-dependent Ginzburg-Landau theory within a Gaussian approximation. 



APPENDIX B: DC FLUX TRANSFORMER EXPERIMENTS 



DC flux transformer experiments have been used as probes of the nonlocal resistive response of iiortex liquids. In 
particular, by comparing transport measurements in twinned and untwinncd YBCO, Lopez et al£3 argued that in 
the twinned crystal, where the liquid transforms into a solid via a continuous glass transition, measurements in the 
transformer configuration may provide a direct probe of the longitudinal correlation length near the transition. Here 
we reconsider the analysis of such experiments in the context of hydrodynamics following earlier work by Huse and 
Majumdar.E3. Our objective is to make the connection between the measured voltages and the diverging correlation 
lengths more precise. As pointed out by Huse and Majumdar, explicit calculations of the current and voltage patterns 
have to be carried out for each specific sample geometry by solving a fourth order partial differential equation with 
rather complicated boundary conditions. Here we reformulate the problem discussed in Ref. |3^ incorporating all 
relevant viscosity coefficients, then discuss the solution in terms of diverging length scales. We also correct Eq. (7) of 
Ref. ^ where some boundary terms were left out. 

We are interested in discussing dc transport experiments in finite materials of specified geometry, where the current 
is injected and withdrawn via contacts placed at the materials' boundaries. As discussed by Huse and Majumdar, 
in the bulk of the sample in steady state the total current satisfies V • j = and the electric field is irrotational, 
V X e = 0. The relationship between cu rren t and field is given by the nonlocal Ohm's law, Eq. (2.5), which at long 



wavelength takes the form given in Eq. (Al). The field e can be written in terms of a local voltage as 



(Bl) 



and the nonlocal Ohm's law transforms into a fourth order differential equation for the voltage, $, 



(^p,^iO)dpd^ - Sf,ap^df,dadpd^ $ = 



(B2) 
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The sample boundaries separate a region where the local field or vortex flow velocity can have spatial gradients 
due to viscous forces among the vortices from the outside region where no such spatial gradients occur. Neglecting 
nonlocalities on length scales of the order of the penetration length, the first spatial derivatives of the field e (or 
equivalently of the vortex flux tensor Qfj,^) can therefore exhibit jump discontinuities at the sample boundaries. 
Such discontinuities correspond to (5-function contributions to the local current, denoted by j"^. Considering a local 
coordinate system with an s axis normal to the sample boundary, located at s = 0, and directed inwards, the nonlocal 
Ohm's law at the boundary takes the form. 



(B3) 



where the subscript s is not summed over. The last term on the right hand side of the equation equals the S function 
contribution to the surface current. Such contribution must be flowing within the surface of the sample, i.e., its 
component normal to the boundary must vanish. 



= -Ss 



= 0. 



(B4) 



This provides the first set of boundary conditions for Eq. (B2) The second boundary condition is obtained by requiring 
that the 5- function part of the divergence of the current at the surface equals the total current Is{s) injected at that 
point of the surface. 



fj,± ssiy^fj,± 



+0+ 



= Is{s), 



(B5) 



where again the index s is not summed over and only runs over the two components perpendicular to the surface. 
This requires that the voltage at the surface satisfies. 



s=0+ 



Is{s). 



(B6) 



The last term in square brackets on the left hand side of Eq. ( JBq ) was absent in Ref. [34 

The specific transformer configurations of interest are shown in Fig. 9. We neglect spatial inhomogeneities in the x 
direction (parallel to the flow) and reduce the problem to a two-dimensional one. The sample has finite extent W in 
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the y direction and thickness L in the external field direction. 





Current is injected and withdrawn at the contacts shown placed on the upper and lower surface of the sample. 
Therefore no current is injected at the boundaries normal to the y direction. In order to discuss the role of the various 
viscosity coefficients in controlling the nonlocal response of the vortex array, it is useful to return to our hydrodynamic 
equations in terms of vortex densities and fluxes. In the geometry of interest there arc two nonvanishing components 
of the electric field from flux motion, related to the corresponding components of the vortex flux tensor as 

^yiV^z) = v^{y,z), 

e^iy,z)^^V{y,z). (B7) 



FromEq. ( |2.15| ) in a steady state vortex current is conserved, which requires 

d,V = dyv,, (B8) 

This condition is equivalent to the requirement that V x e = and implies that the components of the vortex flux 
tensor can be derived by a potential, <i>^, according to 
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Vx{y,z) = dy^y, 



where is simply proportional to the voltage defined in Eq. (Bl), 
governing vortex flow in the transformer geometry shown in Fig. 9 are 



no 00 



1\\V = [ifdl + rildl)v + r, + xdyd^v^ ~ ^dyU + 



(B9) 

The hydrodynamic equations 
(BIO) 

J.. (Bll) 



By differentiating Eq. (BIO) with respect to y and Eq. (Bll) with respect to z and making use again of the fact that 
the total current is divergence-free in bulk, we immediately obtain the equation for the potential <!>„, 



= 0. 



(B12) 



This fourth order differential equation has to be solved with the boundary conditions, 



= 0, 



y=o,w 



=0,L 



-f±dy<t>y - r/9^$„ - {r], + 2r/x + Tf)dldy^y 

7iia,$„ - 7,*af - (77, + 277, + 77*)a2a,$„ 



v=o,w 



=0,L 



±^/.(y,z = ±^). 
71000 2 



(B13) 



To proceed, we follow Huse and Majumdar and neglect viscous drag associated with 77 and rjl (the bulk viscosity 
77b drops out of the problem for linear incompressible flow). This amounts to dropping the terms ~ dy and ~ 9^ 
compared to the fourth order cross derivatives in Eq. ( B12| ). Letting 77^ — rjz + rf + 2772;E2l, the equation simply reduces 
to 



iLd'y^v + 711 5,'$. - fl.dpi^y = 0, 



(B14) 



with the boundary conditions 



"f^dy^y - fj^dldy^y 



y=Q.W 



z=0,L 



0, 



71000 



/,(y,±L/2). 



(B15) 



The solution for an arbitrary distribution of current sources and sinks has been given by Huse and Majumdar.EHl 
Simply translating their results into our notation, we write the current sources at the top and bottom of the sample 
as 



/(y, z = ±L/2) = 5^ 7± cos(7i7ry/W^). 

n=0 

The integral of the current injected or withdrawn over the entire sample surface must vanish. This requires 
Making use of the boundary conditions at a; = 0, W , the voltage is then written as 

00 

$(y, z) = ^ Un{z) cos{mTy/W), 

where 



(B16) 



(B17) 



n=0 



Uq{z) = Co + Soz, 

Un{z) — Cn COSh(iir„z) + Sn sinh(i4r„z), 71 > 0, 



(B18) 
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with 



711 (r) 



(mr/Wy 



-1-1/2 



(B19) 



The coefficient Cq sets an arbitrary zero for the voltage. We will let Cq — below. The other coefficients are given by 



So = Pii^o ' 



2 cosh(i^„L/2)' 



n > 0, 
n > 0. 



(B20) 



2 sinh(ii:„L/2) 

Now we use this solution to discuss the experimental finding of Lopez et alS in twinned and untwinned samples 
of YBCO. These authors have carried out two sets of measurements. The first consists in injecting current at the 
left contact at the top of the sample and extracting it at the right contact, as shown in Fig. 9. The voltages Vtop 
and Vbottom across the two central contacts at the top and bottom of the sample are then measured as functions of 
temperature. For this choice of current sources, — 0, for all n, and Iq = 0. The voltage at the top and bottom of 
the sample is then given by 



a>(y, ±L/2) = -T^P^Y. cos(n^2//M^)if„ (— ^ 
2 ^ — ' VriTT 



cosh(i^„L/2) _^ sinh(i^„L/2) 
sinh(ii:„L/2) cosh{KnL/2) 



(B21) 



The temperature dependence of the voltages is controlled by that of the inverse length Kn, given in Eq. ( B19 ), 
which in turn depends on the friction and viscosity coefficients which diverge near the vortex glass and the Bose glass 
transitions with the critical exponents discussed earlier. In particular we recall that 



t 



711 m ll^iw 
711 (T) j^^Kao) ^ ' 



2<.(C-1) 



(B22) 
(B23) 



where we have extracted the divergence by rewriting all transport properties in terms of their bare values (denoted by 
a superscript """), which depend only weakly on temperature at the transition. The longitudinal correlation length 
will be defined as 



We begin by discussing the case of vanishing viscosity. Then A'„ is simply given by 



mr /7_L 



nn / p\\ 



2f(C-l) 



(B24) 



(B25) 



It vanishes therefore at the Bose glass transition {( = 2), but it remains finite at the isotropic vortex glass transition 
= 1). Choosing for simplicity a current distribution with i'^ ^ and i+ = for all other values of n, the voltage 
for fjz — is given by 



TT 



cosh{K^L/2) smh{KfL/2) 



± 



smh{K°L/2) cosh{K^L/2) 



(B26) 



with Ki given by Eq. (B25) for n — 1. Close to the glass transition where Ki ^ t'^'^^'' -'^^ we can use K^L << 1 to 
obtain. 



$0(2/, ±L/2) « -P^^it cos{Try/W) 



(B27) 



and the voltage at the top and bottom of the sample are equal. Conversely, if K^L » 1, we find 
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w . 

^o{y,+L/2) Ri -^/pjjpX— i+cos(7rj//W), 

$o(y,-i/2)«0, (B28) 

with corrections to the above expressions decreasing exponentially as e~^°^. In this case the vortices behave like 

uncorrelated pancakes and the bottom voltage is much smaller than the one at the top. In other words, there is a 
crossover at K^L ~ 1, corresponding to p^^L"^ /A ~ piW^ /tt"^, from a regime where the top voltage exceeds the bottom 
voltage to one where both voltages are the same. Near the Bose glass transition, where t^" , this crossover 

is sharp and occurs at a temperature that depends on both the width W and the thickness Lof the sample and is 

given by t"^ ^ (2W/ ttL) ^^7" / 7^^ . In contrast, near the vortex glass transition where C = !> is weakly temperature 

dependent and no sharp crossover is expected. 

Consider now the response in the case where the viscosity fiz is very large, with (%/7||)(7r/W)^ >> 1, so that for 
all values of n, 

Kn^K^ = .l^=^. (B29) 

V ^11 

In this case the top and bottom voltages for the first experiment (Fig. 9) are given by 

W 



$^(2/,±L/2)«-i^ 



cosh(L/2C||) ^ sinh(L/2^||) 



^i+cos(n7ry/W)(— ) . (B30) 



n=l 



_sinh(L/2^l|) cosh(L/2^||) 
The crossover is now controlled by the length ^|| t"^". If i/2^|| >> 1, the voltage reduces to 

$^(y,+L/2)«-^ Vz+cos(n7ry/VF)(— ) , 

$oo(2/,-i/2)~0, (B31) 

i.e., voltages are uncorrelated well above the glass transition. As the glass transition is approached and i/2^|| << 1, 
we obtain, 

$00 (y, ±L/2) ^-E±J2it cosinwy/W) (—) . (B32) 

n— 1 

In this case the top and bottom voltages are the same and they are both reduced by a factor ^\\/L compared to the top 
voltage at high temperature. Their vanishing at the transition is controlled by the vanishing of p± . We find a sharp 
crossover from a high temperature region where the top voltage exceeds the bottom voltage to a low temperature 
region where the two voltages are equal, the crossover temperature T* depends on the sample thickness and is given 
by ^\\{T*) ~ L/2. It is therefore expected to decrease with sample thickness according to T* ^ TBG{L/aQ)~^/^" , 
consistent with the observations of Lopez et al. in twinned samples. For T < T* the voltages are the same and 
decreasing with decreasing temperature, vanishing at Tbg as p±^ ~ t"'^''~'^\ 

The second set of measurements carried out by Lopez ct al. consists of injecting a current at the top surface and 
extracting it at the bottom surface, while measuring the voltage across two contacts &t z = ±L/2 and at the same y 
position, as sketched in Fig. 9 . The corresponding distribution of current sources and sinks corresponds to z+ — 
for every n. The voltage difference between top and bottom is given by 

A$(y) = $(2/,L/2)-$(y,-L/2) 

= iQP\\L + 2p^\^i-cos{n'Ky/W)Kn\^ — j tanh(/ir„L/2). (B33) 

For simplicity we assume the currents injected at the top and withdrawn at the bottom are uniform in y, i.e., 
Iq = —Iq = constant, and = — i+ = for n 7^ 1. the voltage is then simply 

A<I)(y) = Zo-p||L (B34) 

and for all values of y it is controlled by the c-axis resistivity py, which vanishes at the glass transition. 
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